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If two identical emitters are coupled to a common reservoir, entanglement can be generated during the 
decay process. When using Bell's inequality to examine the non-locality, however, it is possible that the 
bound cannot be violated in some cases. Here, we propose to use the steering inequality to examine 
the non-locality induced by a common reservoir. Compared with the Bell inequality, we find that the 
steering inequality has a better tolerance for examining non-locality. In view of the dynamic nature of the 
entangling process, we also propose to observe the quantum coherent dynamics by using the Leggett-Garg 
inequalities. We also suggest a feasible scheme, which consists of two quantum dots coupled to nanowire 
surface plasmons, for possible experimental realization. 



Spontaneous emission (SE) occurs when a two-level emitter exponentially relaxes into its ground state and is 
one of the fundamental concepts in quantum mechanics. Its history can be traced back to such early works 
as that of Albert Einstein. When considering an ensemble of two-level emitters, the emission rate can be 
greatly enhanced, a phenomenon known as superradiance (SR) 1 . The phenomenon of SR has been widely studied 
in atomic 2 and solid state 3 " 5 systems. Its role and applications in quantum information science have also been 
investigated recently 6 ' 7 . The key ingredient in SR is that the emitters exchange excitations with one another via the 
common reservoir. On one hand, this process creates non-local quantum coherence (entanglement) between the 
emitters 8 . On the other hand, one should also note that it is a time- dependent process which involves quantum 
coherent dynamics. 

To test the non-locality, the Bell inequality is probably the mostly utilized tool 9 . More recently, an inequality 
based on the concept of quantum steering was proposed 10 . This enables one to make a comparison between the 
hierarchy of states which violate Bell's inequality, quantum steering, and entanglement. Focusing on a single 
object, Leggett and Garg in 1985 derived an inequality 12 to test the assumptions of macroscopic realism and the 
possibility of non-invasive measurement. An experimental violation of this inequality in a superconducting 
circuit has recently been seen 13 . The Leggett-Garg (LG) inequality can also be applied to microscopic sys- 
tems 14 " 18 to examine the quantum coherent dynamics therein. 

The entanglement generated during SR (or reservoir-induced entanglement) is a well-known phenomenon. 
However, to our knowledge, there are no studies on exactly how well the states generated by SR, in a realistic 
situation, can violate the hierarchy of inequalities outlined above. Actually, the superradiant process provides a 
good platform to examine not only non-locality, but also the quantum coherent dynamics. Therefore, we will 
separately investigate these two quantum features of the simplest superradiant case: SE of two identical emitters 
coupled to a common environment. As is well known, and will be demonstrated explicitly here, the entanglement, 
the steering inequality, and the Bell inequality form a logical hierarchy. Since the collective decay is a dynamical 
process, we then propose to use the LG inequalities 12141519 to examine the quantum coherent dynamics. We will 
show that the LG inequality can be violated when choosing a suitable observable. We will also utilize the recently- 
developed extension of Leggett-Garg inequality (ELG) 14,15,19 to further study the quantum coherent dynamics. For 
the experimental realization, we will provide a feasible scheme, which consists of two quantum dots coupled to 
nanowire surface plasmons 20 . 
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Results 

Consider two identical two -level emitters coupled to a reservoir. The 
interaction Hamiltonian is written as 



H int =Yj H [-S (*! U k + °2 a k eitd ) + H ' C - 



(1) 



where a~ k is the annihilation operator of the A:- mode field in the 
reservoir, is the raising operator of the jth two-level emitter, 

and d is the inter-emitter distance vector. Here, g% is the coupling 

strength between the jth emitter and the /c-mode field. 

In the interaction picture, the evolution of the total system 
(denoted by the density matrix x) is governed by the Schrodinger 

equation % = — [H int ,%] . After performing the Born-Markov approxi- 
ih 

mation, the equation of motion of the reduced density matrix p in the 
single excitation manifold reads 
y Sm(qd) 
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where d = 1 5 , p(t) = Tr R [#(£)] with R denoting the reservoir, and y is 

the spontaneous emission rate of the emitter. Here, q = co 0 /c with co 0 
being the two -level energy spacing of the emitter and c being the 
speed of light. 

In the single excitation manifold, p(t) can be spanned by the basis 
{| 1, 0), |0, 1), |0, 0)}, where 1 1, 0) (|0, 1» denotes the first emitter is in 
the excited (ground) state with the second one in its ground (excited) 
state, and |0, 0) represents both the emitters are in the ground state. 
The inset in Fig. 1 shows the quantum coherence (|1, 0) (0, 1|) 
between the two emitters in the limit of d = 0. As seen, the coherence 
saturates to a fixed value. This quantum coherence, created by the 
cross terms in the master equation, leads to the entanglement 
between the two emitters. 

Quantum entanglement induced by the common reservoir. To 

demonstrate the degree of the entanglement, we use p to calculate 
the concurrence 21 C = C(p s ) = maxf^ — X 2 ~ ^3 ~~ where A t -s are 
the square roots of the eigenvalues of pp with p= ((r y ®(Ty) 
p* (o y ®(Jy). In Fig. 1, we assume two different initial states: |1, 0) 

for Fig. 1(a) and the triplet state (|1,0) + |0,1» for Fig. 1(b). As 

V 2 

shown by the red-dashed curve in Fig. 1(a), the excitation is 
transferred to the second emitter with the coherence induced by 
the common reservoir. The entanglement is created during the 
excitation transfer and becomes saturated as t —> 00 . However, 
when the initial state is the triplet one, the entanglement (red- 
dashed curve) decays exponentially with an enhanced rate as 
shown in Fig. 1(b). 

It is interesting to compare the concurrence to the quantum dis- 
cord 22 , which is a measure of quantum correlations between two 
subsystems of a quantum system. In our system, if we denote the 
first (second) emitter as subsystem A (B), the quantum mutual 
information can be defined as 

I(A:B) = S(p A )+S(p B )-S(p), (3) 

and its equivalent expression J [A : B) can be written as 

J m (A:B) = S(p A )-S m (A\B). (4) 

The quantum discord V(A : B) is then defined as 22 ' 23 

V(A : B) =X(A : B)-J(A : £), (5) 
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Figure 1 | The concurrence and quantum discord. The concurrence (red- 
dashed curve) and the quantum discord (black-solid curve) of two-emitter 
superradiance for the initial states: (a) 11, 0) and (b) —= (|1,0) + 10,1)). The 

complete set of the projectors {IT;}{cos0l 1) + e^sin^lO), e ^sin0l 1) — 
cos^l 0)} is found to contribute maximally to the value of t7{n,-} (A '• f° r 
both initial states when 6= — (for any (/)). Inset: The coherence (off- 
diagonal) term, 1 1, 0) (0, 1 1. In plotting this figure, we have assumed the 
inter-emitter distance: d = 0. 

where S(p) = — Tr(plogp) (in what follows, we use the natural 
logarithm) is the von Neumann entropy of the system, Pa(Pb) is 
the reduced density matrix of p after tracing out B (A), and 
S{u.}(A\B) in Eq. (4) is the quantum conditional entropy describing 
the average uncertainty about the state of A given that the state of B is 
known. J^{A : B) is generalized from its classical counterpart 
through measuring the subsystem B by a complete set of projectors 
{IT f }. Since J{Ui} (A : B) strongly depends on the projectors {IT,-}, the 
quantity J {A : B) in Eq. (5) denotes the maximum of J^ Ui }(A : B) 
over all {II,}. 

We plot the quantum discord by projecting the subsystem B with 
the complete set {cos#|l) + £^sin#|0), e -i *sin0|l) - cos#|0)} for the 

initial states |1, 0) [Fig. 1(a)] and (|1,0) + 10,1» [Fig. 1(b)], 

v2 
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respectively. As shown by the black-solid curves in Fig. 1, though the 
quantum discord has been proven to be more robust 24 than the 
concurrence in some systems, the behavior of the quantum discord 
resembles that of concurrence. This is because the entanglement is 
always present in our system, and thus there is no such thing as, 
"quantum discord without entanglement" here. In this case, the 
concurrence and the quantum discord reveal similar behavior. 

Testing the non-locality. Let us first utilize the Bell inequality to test 
the non-local properties of the emitters during the collective decay. 
The Bell quantity associated with the CHSH inequality has the 
following form 25 , 

BcHSH = a-(j® (h + b'^'G + a'-d® (b — b^j-G, (6) 

where a, a', b, V are unit vectors in R 3 . Here, a-G= YM=\ a i G ^ wnere 
Gi is the standard Pauli matrix. Then, the CHSH inequality of a state 
p is 

\(BcHSH) p \ = \Tr(pBcHS H )\<2. (7) 

The maximum value of the CHSH inequality is given by 

B m ax= max Tr(p5cffSff). (8) 

a,a',b,b 

Assuming a general initial state: cos a 1 1, 0) — sin a|0, 1), we plot in 
Fig. 2(a) the maximum value of the CHSH inequality (B max ) as 




Figure 2 | Testing the Bell-CHSH and steering inequalities during the 
superradiant process. Assuming a general initial state: cos a 1 1, 0) — 
sin a 1 0, 1), we plot in (a) the maximum value of the Bell-CHSH inequality 
and in (b) the steering parameter as functions of time and a in the limit of 
d = 0. In plotting the figure, the time t is in unit of 1/y. 



functions of time and a in the limit of d = 0 for the superradiant 
process. 

As a comparison, let us also consider the steering inequality. If the 
correlation measurement on the emitter- 1 (B) and emitter-2 (A) is 
described by the probability distribution P[B t = b, A t = a) with b = 
± 1, and a = ± 1, the steering inequality is written as 26 

N r 2i 

s^^eUa^Ui, (9) 

i=l 

where N(= 2 or 3) is the number of mutually unbiased measure- 
ments (the Pauli X, Y and Z) for the emitter-2, and 

b=±\ 

with emitter-2 's expectation value for a measured (conditioned on 
emitter- l's result) defined as 

(Ai) Bi = b = P(A t = + 1 \Bt = b)- P(A t =-\\Bi = b). (11) 

In Fig. 2(b), we also plot the steering parameter (S 3 ) as functions of 
time and a in the limit of d = 0. 

As shown in Fig. 2, we can see that the area of violating the steering 
inequality is larger than that of violating the Bell inequality. In the 
limit of a —> 0 (only one of the emitters is excited), Bell's inequality 
can never be violated, while the steering one can still be violated. This 
again shows that the steering inequality has better tolerance in exam- 
ining the non-locality of the entangled states 10 , as expected. Together 
with the results in Fig. 1, one can see that the entanglement, the 
steering inequality, and Bell's inequality form a logical hierarchy 11 , 
i.e. every violation in Bell's inequality leads to the violation of the 
steering inequality (but not vice versa) and every violation in the 
steering inequality is a entanglement witness (but not vice versa). 
The underlying physics of this logic hierarchy can be understood as 
follows. From the operational definition 11 of Bell non-locality and 
steerability, Bell non-locality can be viewed as a three-party task: 
Alice and Bob, who cannot communicate with each other, try to 
convince Charlie that the state they share is entangled. Charlie does 
not trust them and can be convinced only if the correlations between 
them do not fit the local hidden variable model. However, the steer- 
ability is regarded as a two-party task: Alice tries to convince Bob that 
the state she prepares is entangled, but Bob does not trust her. Bob 
can be convinced only if the correlations between his local measure- 
ments and Alice's results cannot be explained by a local hidden state 
model. In Bell's non-local scenario, if we relax the condition to make 
Bob trustworthy, the result will be consistent with that of the steering 
one. Bell's non-locality is therefore a stronger non-local concept 
compared with the steerability. In other words, an entangled state, 
which violates the steering inequality might not violate Bell's 
inequality. The steering inequality therefore has better tolerance in 
examining the non-locality. 

Testing the quantum coherent dynamics. Given an observable Q(t), 
which is bound by \Q(t)\ < 1, the assumption 12 of realism and 
noninvasive measurability implies the Leggett-Garg inequality 12 : 

LG Q = (Q(fi)Q) + (Q(h + fe)Q(fi)) - (Q(h + h)Q) < 1, (12) 

here Q = Q(t = 0), and t x < t 2 . The experimental violations of this 
inequality in a "macroscopic" superconducting circuit 13 , photon 
polarization states 27 " 29 , electron-nuclear spin pairs 30,31 have recently 
been seen. 

Normally, the LG inequality is applied to a macroscopic object to 
verify the property of macroscopic realism. Here, we wish to make 
use of it as a tool to check the quantumness (quantum coherent 
dynamics) during the superradiant process. By choosing ti = t 2 = 
t with the initial state being in |1, 0), we apply the original LG 
inequality in Eq. (12) to examine the quantumness in the two-emitter 
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Figure 3 | The original LG inequality for as a function of time. The blue- 
dotted, dashed, and solid curves represent the results of the observables 
Q = Qi= ll,0)(l,0l,Q = Q 2 = 1 1, 0> <1, 01 + 10, 0) (0,01, and Q= Q 3 = 
1-2 10, 1) (0, 11, respectively. In plotting this figure, the initial state is 
assumed to be 11, 0). 

SR scenario. If we choose Q = Q x = |1, 0) (1, 0| as our observable, 
there is no violation of LG Q (t) as shown by the blue-dotted curve in 
Fig. 3. This coincides with our intuition: the decay of the state 1 1, 0) is 
a monotonically decreasing function. However, if we choose Q = Q 2 
= 1 1, 0> (1, 0| + |0, 0) (0, 0|, i.e. the excitation is not in the second 
emitter, the red-dashed curve shows the violation of LGq 2 (t) in the 
early stage of the time domain. We also plot in Fig. 3 the result of Q = 
Q 3 = 1-2 1 0, 1) (0, 1 1 . One can see that the maximum of the violation 
LGq 3 (^) is enhanced. The analytical results of L Q (t) can also be 
worked out. For instance, the expression of LGq 3 (t) is 

LG Q3 (t) = ^+le-^e-^[-2e^e^-l+e^) 2 + (-l + e 2 ^ 2 ],(13) 



where = 



-.l+-e 

Sin(qd) 
qd 



yt. 



In practice, the non-invasive measurements are not easily imple- 
mented. Recently, an extension of the LG inequality was developed 
for an open transport system 1415 . Based on the Chapman- 
Kolmogorov equation in stochastic theory 14 , when measuring the 
population of a state of a classical Markovian system, the following 
inequality can be derived: 



\L Q (t)\ = \2(Q(t)Q)-(Q(2t)Q)\<(Q), 



(14) 



where (Q) is the expectation value of the zero-time population Q = 
Q(t = 0), and (Q(f)Q) is the two-time correlation function. Note that 
if the zero -time state is the steady state then this is equivalent to the 
original 12 LG inequality when the measurements are performed non- 
invasively. Violations of Eq. (14) mean that the dynamics is beyond 
the classical Markovian description. Note that in deriving Eq. (14) 
the assumption of the non-invasive measurement is not required in 
that one can replace the initial measurement by state preparation. 

The blue-dotted, red-dashed, and solid curves in Fig. 4(a) repres- 
ent the results of \L Q (t) | for the observables Q = Q l = | 1, 0) (1, 0|, Q 
= Q 2 = \h 0) (1, 0| + |0, 0) (0, 0|, and Q = Q 3 = 1-2 |0, 1) (0, 1|, 
respectively. Again, we see that the violation depends crucially on the 
choice of the observable. Figure 4(b) is the 3D plot of |Lq 3 (£)| as 
functions of time t and inter- emitter distance qd. As expected, the 
maximum of the violation gradually decays to unity when increasing 
d. As qd = nn, where n is an integer, the function |Z,Q 3 (f)| is not 
violated because the two emitters act as independent objects due to 
the interference effect in SR 32 ' 33 . This can be seen more clear from the 
analytic expression of Lq 3 (t). 
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Figure 4 | Violations of ELG inequality during the spontaneous decay of 
two emitters. With the initial state is initially in 1 1, 0), we plot (a) The ELG 
inequalities for \L Q {t) l/(Q) as a function of time. The blue-dotted, red- 
dashed, and solid curves represent the results of the observables Q = Qi = 
1 1,0) (1,01, Q= Q 2 = 1 1,0) (1,01 + 10, 0) (0,01, and Q= Q 3 = 1-2 10, 1) 
(0, 11, respectively, (b) The ELG inequalities for |Lq 3 (f) | / (Q3) as functions 
of time and inter-emitter distance, (c) The border line (|Lq 3 | = 1) that 
separates the violation and non-violation regime. As qd = nn, \Lq 3 | is equal 
to unity as shown by the red-dashed lines. 

As pointed out in Ref. 14,15, the ELG inequality is equal to the 
original LG inequality if (i) the zero -time state is the steady state, (ii) 
the operator Q in Eq. (9) is normalized as Q = 2Q/(Q) — 1, and (hi) 
the measurements are performed non-invasively. For the observable 
Q 3 , surprisingly, we find that the kernel is mathematically identical to 
that of the LG inequality, i.e., 

L Q3 (t)/(Qs)=LG Q3 (t)- (15) 

However, in plotting Fig. 3, we have assumed the excitation is ini- 
tially in the state 1 1, 0), not the steady state. To answer this, we check 
the two-time correlation functions and find that (Q(f)Q) = 
<Q(2f)Q(f)> for the observable Q = Q 3 . This explains why the ELG 
inequality is mathematically equivalent to the original LG inequality. 
However, one should note that although the mathematical results are 
the same, it does not mean that they are physically equivalent since 
the originally LG inequality requires the assumption of non-invasive 
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measurement. We can only say that the violation of the ELG inequal- 
ity is the indication of the quantum coherence. 

It is also interesting to see how the time interval of violation (t^J 
varies with inter-emitter distance qd. We therefore plot in Fig. 4(c) 
the border line (|Lq 3 | = 1) that separates the violation and non- viola- 
tion regime. The interval of violation is maximum when d = 0 and 
gradually decreases to a fixed value with the superradiant interfer- 
ence feature. To explain the approach to the fixed value, we take the 
series expansion of Eq. (10) to the order of fi 2 . Lq 3 is then approxi- 
mated as 



excite 



1 



2yt 



(4-2e yt )p 2 . 



(16) 



Therefore, the time interval of violation is T vio = In 21 y in the limit of 
large d. Of course, as d— > o°, the value of |Lq 3 | should go to unity since 
P — > 0 in this limit, i.e. no superradiant effect. Equation (1 1) is to tell 
us that T vio is kept fixed, not like the maximum of the violation, when 
the inter- emitter distance d gradually approaches to infinity. 

A proper choice of the observable is the key to see the violation of 
the ELG inequality. This can be understood by the following 
example. Consider a qubit in a lossy cavity. Let us choose | + , 0 C ) 
( + , 0 C | or | — , 1 c ) ( — , 1 c \ as our observables, for which we plot in Fig. 5 
the function \L Q \. Here, the state | + , 0 C ) (| — , l c )) denotes the qubit is 
in its excited (ground) state and there is zero (one) photon in the 
cavity. For the observable | — , l c ) (— , l c \, it is possible that \L Q \ is not 
violated (the red-dashed curve) even though the population 
dynamics is still oscillatory. Therefore, it is very important to choose 
a proper observable when using the ELG inequality as a tool to 
indicate the quantumness. 



Discussion 

Here, we propose an experimental realization for verifying the 
quantum coherent dynamics in SR. Consider two quantum dots 
positioned near a metal nanowire 20 ' 34 ' 35 as shown in Fig. 6. Due to 
the quantum confinement, the surface plasmons propagate one- 
dimensionally along the axis direction on the surface of the nanowire. 
Quantum dot-1 is assumed to be excited initially. In one-dimen- 
sional problems, the SR becomes independent of the inter-dot 
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Figure 5 | The ELG inequality for Rabi-oscillations between a qubit and 
cavity photon. With the state I + , 0 C ) being initially excited, the red curves 
are the results of the observable I — , 1 c ) {— , 1 J , while the black curve is the 
result of the observable I +, 0 C ) ( + , 0 C I . Inset: Population dynamics of the 
state I — , 1 c ). Here, k and y are the cavity loss and atomic dissipation rate, 
respectively, while g is the coupling strength between the qubit and the 
cavity field. 



QD-1 
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Figure 6 | Schematics of a experimentally-accessible system. Two two- 
level quantum dots coupled to metal nanowire surface plasmons. Here, the 
quantum dot- 1 is assumed to be excited initially, and one measures the 
populations of the quantum dot-2. 



distance. If we choose Q = Q 3 , the kernel of both LG and ELG 
inequality can be written as 

L Q3 (t)/{Q 3 )=LG Q3 (t)=\(l + e- i " t -4e- 2 y' + 4e-y t ). (17) 

The advantage of choosing Q = Q 3 = 1-2 |0, 1) (0, 1| is that, experi- 
mentally, one only needs to measure the population of quantum 
dot-2. 

Summarizing, we have utilized both the Bell inequality and the 
steering inequality to examine the non-locality during the superra- 
diant decay and have shown that the steering inequality has better 
tolerance in examining the non -locality. Moreover, we have also 
pointed out that the quantum coherent dynamics during the super- 
radiant decay can be examined by using the LG inequality. Although 
violating one of the inequalities indicates that quantum effects exist 
during the superradiant process, the violations have two different 
meanings: one is for the non-locality (Bell, steering), and the other 
is for the quantum coherent dynamics (LG). In experiments, we 
stress that only the measurements on one of the emitters are required 
to test the LG inequalities. The ELG and LG inequalities are math- 
ematically equivalent when choosing a proper observable. Therefore, 
the ELG inequality, which does not require the non-invasive measur- 
ability, can in practice assist in verifying the quantumness in SR. As a 
final remark, we wish to point out that the contextuality test 36 can 
also be used to indicate the non-locality mentioned above. Similar to 
the LG inequality, its temporal scenario 37 can be utilized to indicate 
the quantum coherent dynamics. Further investigations in this dir- 
ection are under way. 

Methods 

Details of the derivation of Bell-CHSH inequality. To obtain the results in Fig. 2(a), 
we adopt the method provided in Ref. 26. In order to verify whether the Bell-CHSH 
inequality is violated for a given state p, one has to maximize the mean value of the 
Bell-CHSH parameter: 

(18) 



B max = max Tr(pBcHS H ) =2y/M(p) = 2y% +^2, 

a,a',b,b 



where ki and k 2 are the two largest eigenvalues of TjT p . Here, T T p is the transpose of 
T pi and the matrix elements of T p is written as 

t mn =T:r(pG m ®G n ), (19) 

where o n is the standard Pauli matrix. 

Details of the derivation of steering inequality. To obtain the results in Fig. 2(b), we 
first measure the state p in the z-basis, and the probability that Bob obtains ± 1 is 

p(Bi= + l) = Tr(nf i = ±1 p), (20) 



where 



and 



1 0 
0 1 



1 0 
0 1 



<8> 



1 0 
0 0 



0 0 
0 1 



(21) 



(22) 



After Bob's measurement, the state is mapped to 
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Tr nf 



With the state p' , the probability that Alice obtains ± 1 is 
P(A i =±l\b)=Tr(fl*' = ±1 p'). 



(23) 



(24) 



With Pfa = ± 1) and P{A t = ±l\b), Eqs. (10) and (11) are obtained accordingly. To 
calculate S 3 in Eq. (9), one repeats the above procedure by replacing Ylf = ±l , ft?' = ± 1 
with fl J = ± \tlf = ±l and ft^ = 1 1 , tif = ± 1 . One can then have S 3 in three mutually 
unbiased basis. 
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